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RSRE  MEMORANDUM  NO  3355 


MOMENT  METHOD  ANALYSIS  OF  VHF  ANTENNAS  ON  VEHICLES  ON  AN  IMPERFECT  GROUND 
D  Williams  &  D  J  Brammer 

INTRODUCTION 

'  The  study  of  vhf  vehicle  antennas  Is  complicated  by  the  Irregular  vehicle 
which  Is  typically  of  the  order  of  a  wavelength  and  the  Imperfectly  conducting 
ground  on  which  the  vehicle  stands.  This  memo  presents  a  treatment  of 
Idealized  armoured  vehicles  fitted  with  whip  antennas  standing  on  an 
Imperfectly  conducting  ground  based  on  the' Moment  Method' for  the  vehicle  and 
Including  a  Sommerfeld  based  treatment  of  the  grotmd.  The  20-76  MHz  vhf  band 
Is  of  Interest  but  since  the  Investigation  showed  that  ground  effects  are  more 
Important  at  the  lower  part  of  the  band  attention  was  focussed  on  this. 

■.  r  -  I  ’  ' 

Idealized ^versions  of  the  armoured  vehicles  of  Interest  were  examlned|> 

Fig  1  shows  the  simplest  -  a  rectangular  box;  Fig  5  shows  a  shape  nearer  the 
true  vehicle.  The  Idealized  models  are  fitted  with  2M  whip  antennas  In  the 
position  used  In  practice.  The  vehicles  are  8M  long  and  are  supposed  to  be 
supported  by  Insulating  tracks  30cm  above  ground.  The  numerical  analysis  was 
supplemented  by  some  experimental  work  carried  out  at  1/30  scale  In  an  anecholc 
chamber. 

MOMENT  METHODS 

There  are  many  variations  of  the  Moment  Method  ranging  from  the  basic 
formulation  to  details  of  the  numerical  analysis.  A  wire  grid  modelling  of  the 
surface  of  the  vehicle  Is  the  most  common  approach.  The  wire  model  of  Fig  1, 
which  Includes  280  segments  was  the  subject  of  the  first  Investigation.  The 
methods  tried  were: 

1  Pulse  current,  point  matching  (collocation). 

11  Piecewise  Sinusoidal  with  Galerkln  Matching  (PWSG)  (1) . 

Ill  A  3  term  (dc,  sinusoidal,  co-slnusoldal)  current  expansion  on 

segments  with  current  and  slope  of  current  continuity  at  segment 
junctions  as  used  by  Burke  and  Pogglo  (2). 

Iv  A  simple  surface  patch  technique  (2). 

In  the  first  3  (wlregrld)  methods  the  Electric  Field  Integral  Equation  Is 
used,  this  Implies  the  use  of  the  tangential  electric  field  In  the  boundary 
condition  whereas  In  the  surface  patch  technique  the  Magnetic  Field  Integral 
Equation  Is  used  with  the  magnetic  field  boundary  condition.  In  the  wlregrld 
analysis  a  difficulty  occurs  with  the  diameter  of  wire  used  In  the  grids 
modelling  surfaces.  The  diameter  chosen  was  In  terms  of  equivalent  surface 
area  selected  by  comparison  of  experimental  results  and  calculations  of  the 
polar  patterns  of  spheres  with  whip  antennas.  The  first  2  programs  were 
written  locally.  The  second  pair  were  available  on  the  program  NEC  (2) 
written  at  the  Lawrence  Livermore  Laboratory  of  the  University  of  California. 
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Fig  2  shows  the  polar  diagram  of  the  vehicle  model  of  Fig  1  at  35  MHz  as 
calculated  by  the  various  moment  method  systems.  The  first  three  methods, 
pulse  current  point  matching,  the  piecewise  sinusoidal  with  Calerkln  matching  and 
the  3  term  NEC,  all  use  288  segments.  The  patch  method  used  patches  plus  wire 
antenna  segments.  Also  shown  is  an  experimental  curve  obtained  in  an  anecholc 
chamber  with  a  continuously  skinned  box  scaled  1/30.  It  can  be  seen  that  all 
methods  give  reasonably  good  results;  the  pulse  current  method  being  slightly 
poorer.  The  piecewise  sinusoidal  method  was  rejected  since  on  modelling  a 
surface  with  a  wiregrid  for  a  given  number  of  segments,  the  piecewise  sinusoidal 
method  Involves  an  Increase  in  the  number  of  unknowns  by  a  half.  The  solution 
of  a  system  of  simultaneous  equations  is  central  to  the  moment  method  and 
increasing  the  number  of  unknowns  by  a  half  Increases  the  niimber  of  operations 
by  a  factor  of  3.4.  The  NEC  3  term  wiregrid  program  was  used  for  the  subsequent 
work.  The  program  contains  many  options  including  the  facility  for  inserting  a 
reflexion  coefficient  approximation  imperfect  earth. 

The  moment  method  solution  gives  as  a  primary  product  the  current 
distribution  on  the  conductor  system  and  therefore  estimates  the  input 
admittance.  Generally  it  is  found  that  the  susceptance  which  is  concerned 
with  near  fields  near  the  feed  point  is  difficult  to  estimate.  Table  1 
contains  estimates  of  input  Impedance  and  is  compared  with  network  analyser 
measurements  on  the  1/30  scale  model. 

TABLE  1 


Freqn 

MHz 

Impedance  Ohms 

Experimental 

Wiregrid 

NEC 

30 

26.1  -  61.5  j 

34.4  -  87.0  j 

35 

48.8  +  3.0  j 

50.5  +  5.6  j 

40 

61.5  +  79.0  j 

71.5  +  87.3  j 

The  results  are  stable  with  respect  to  changes  of  segment  length  as  is 
shown  by  a  series  of  calculations  at  35  MHz. 

TABLE  2 


Segment  Length 

Impedance 

Ohms 

X/20 

50. 5+5. 6 j 

X/28 

53.1+4. 5j 

X/48 

57. 3+2. 6 j 

This  indicates  that  the  current  continuity  algorithm  used  is  reliable. 
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IMPERFECT  GROUND  EFFECTS 


In  mobile  vehicle  to  vehicle  communication  the  vertical  field  at  the 
ground  Is  the  Important  parameter.  The  presence  of  the  ground  has  two  effects 
on  this.  In  the  first  place  the  current  distribution  on  the  vehicle  Is  affected 
by  the  ground.  In  the  second  place  given  a  current  distribution  on  the  vehicle 
the  far  field  Is  affected  by  ground  reflection.  In  the  moment  method  the  field 
due  to  the  current  In  one  segment  Is  to  be  calculated  at  a  match  point  usually 
at  the  mid  point  of  another  segment.  In  the  present  problem  the  distance 
Involved  and  that  of  segment  and  match  point  from  the  earth's  surface  are 
fractions  of  a  wavelength.  This  makes  use  of  Sommerfeld’s  theory  for  fields 
near  the  earth's  surface  necessary.  When  the  observation  point  and  current 
element  are  both  on  the  same  side  of  the  Interface  (In  air)  the  field  Is 
obtained  In  terms  of  the  U  and  V  integrals  together  with  source  and  Image 
terms  (Banos  (3)). 


“  Y-(h+2) 

U  =  f  _  H 


(Xp)X  dX 


'22 


»  -Y^Ch+z) 

/  - 3 — 

'^1^2  ^  ^^2^1 


H  ^^^(Xp)X  dX 
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Where  p  Is  the  separation  of  source  and  observation  points  projected  on  the 
surface  of  the  earth  and  (h+z)  Is  the  sum  of  the  heights  of  these  points  above 
the  earth  (or  the  vertical  separation  of  the  observation  point  from  the  source 
point  Image  in  the  earth's  surface).  The  suffix  1  refers  to  the  ground  and 
2  to  the  air  above,  k  Is  the  propagation  constant. 
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and  H 
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is  the  Hankel  function. 


If  calculations  of  U22  and  V22  were  carried  out  for  each  estimation  of  the 
field  an  Intolerable  computer  processor  time  would  be  taken.  It  was  pointed  out 
by  Brlttlngham  (A)  that  since  with  source  and  observer  points  on  air  the  field 
Is  a  function  of  h+z  and  p  a  two  dimensional  Interpolation  can  be  used.  In  the 
present  work  the  reflection  coefficient  method  of  introducing  the  earth  is  used 
and  the  fields  obtained  are  corrected  by  a  two  dimensional  Interplatlon  from  a 
precomputed  and  stored  table  of  correction  factors  derived  from  the  Lytle  and 
Lager  method  of  applying  the  Sommerfeld  theory  (5). 


Table  3  gives  the  relative  timing  for  the  calculation  of  the  5  field 
components. 
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TABLE  3 


PROCESSOR  TIME 


Free  Space 

Reflection  Coefficient 
Sommer f eld 
Interpolation 


2  msecs 

5.5  msecs 

2.5  secs 

9.5  msecs 


The  wlregrld  model  of  Fig  1  was  analysed  with  the  structure  0.3M  above 
ground  (Insulating  tracks  are  assumed)  with  constitutive  parameters  =  9.0, 

a  "  O.OIS  firstly  at  35  MHz.  The  resulting  polar  diagram  on  the  earth's 
surface  is  shown  in  Fig  3  and  there  is  little  difference  between  the  results 
if  current  distribution  Is  calculated  using  the  S-'^merfeld  treatment,  the 
reflection  coefficient  treatment  or  with  no  earth  at  all.  The  chief  effect 
of  the  earth  Is  noticed  In  calculating  the  polar  diagram  from  the  current 
distribution  and  the  asymptotic  Norton  theory  is  required. 


When  the  same  calculations  are  repeated  at  20  MHz  a  greater  difference  Is 
noted,  the  case  for  the  full  Sommerfeld  treatment  is  more  noticeable  in  Fig  4. 

To  show  that  the  results  do  not  depend  on  the  unrealistic  vehicle  model 
used  and  that  the  current  system  is  not  restricted  to  the  top  of  the  vehicle, 
the  more  realistic  model  of  Fig  5  was  also  analysed.  Again  the  results  in 
Fig  6  show  that  at  35  MHz  the  polar  diagram  at  the  earth's  surface  again 
chiefly  depends  on  the  use  of  the  Norton  asymptotic  theory;  the  method  of 
calculating  the  current  distribution  not  being  Important. 


CONCLUSION 

It  is  concluded  that  the  Sommerfeld  ground  treatment  is  only  important 
at  the  very  lowest  part  of  the  vhf  band  when  calculating  coverage.  When  the 
full  ground  treatment  is  required  substantial  savings  in  computer  processor 
time  are  to  be  found  by  using  the  interpolation  method. 
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FIG.  I  WIREGRID  BOX  REPRESENTATION 
OF  ARMOURED  VEHICLE 
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FIG.  2  RADIATION  POLAR  PATTERN  FOR 
PULSE  POINT  MATCHING,  PWSG.THREE  TERM 
EXPANSION,  PATCH  MFIE  AND  EXPERIMENT 
AT  35  MHz,  (BOX  I  WAVELENGTH  LONG) 
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FIG.  3  POLAR  DIAGRAM  FOR  FREE  SPACE 
CURRENTS  WITH  NO  EARTH  .  FREE  SPACE 
CURRENTS,  RE F^- COEFFICIENT  CURRENTS, 
SOMMERFELD  CURRENTS  ALL  WITH  NORTON 
FAR  FIELD  TREATMENT  35  MHz  (I  X  BOX) 
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FIG.  4  POLAR  DIAGRAM  FOR  FREE  SPACE 
CURRENTS  WITH  NO  EARTH,  FREE  SPACE 
CURRENTS,  REF*-'  COEFFICIENT  CURRENTS, 
SOMMERFELD  CURRENTS  ALL  WITH  NORTON 
FAR  FIELD  TREATMENT  20MHzl8AB0X) 


FIG.  5  SECOND  WIREGRID  REPRESENTATION 
OF  ARMOURED  VEHICLE 
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FIG.  6  POLAR  DIAGRAM  FOR  FREE  SPACE 
CURRENTS  WITH  NO  EARTH.  FREE  SPACE 
CURRENTS,  REfL-  COEFFICIENT  CURRENTS, 


SOMMERFELD  CURRENTS  ALL  WITH  NORTON 
FAR  FIELD  TREATMENT.  35  MHz.  SECOND 
REPRESENTATION  OF  VEHICLE  (I  X  LONG  ) 


